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I.  INTRODUCTION 


In  this  report,  the  formalism  for  calculation  of  the  atmos- 
pheric turbulence -induced  power  spectrum  of  beam  wander  for  a laser 
transmitter  will  be  set  up.  A procedure  for  generating  a random 
sequence  of  values  corresponding  to  each  component  of  beam  wander  will 
be  established.  The  effect  of  atmospheric  turbulence  on  a laser  target 
designator  will  be  evaluated. 

The  temporal  statistics  of  the  effects  of  atmospheric  turbulence 
on  the  direction  of  propagation  of  a laser  beam  is  determined  by  the 
laser  wave  length,  \,  the  laser  transmitter's  (i.e.,  the  laser  target 
designator's)  aperture  diameter,  D,  the  range  from  the  transmitter  to 
the  target,  R,  the  distribution  of  the  optical  strength  of  turbulence 
along  the  path  as  measured  by  the  refractive-index  structure  constant, 

2 

, and  by  the  distribution  along  the  propagation  path  of  the  component 

of  the  effective  wind  velocity  perpendicular  to  the  path,  V ,f.  All 

2 ^ ^ ^ 
of  these  except  CN  and  are  single  value  parameters.  However, 

CN2  and  can  vary  with  the  position  z along  the  propagation  path. 

It  is  to  be  noted  that  the  effective  wind  velocity,  ve££»  is  not  deter- 
mined entirely  by  the  actual  wind  velocity,  but  may  contain  a component 
which  depends  on  the  angular  rate  at  which  the  laser  beam  is  being 
slewed  through  the  air.  What  really  matters  is  the  rate  at  which  a 
"parcel"  of  turbulent  atmosphere  passes  through  the  laser  beam. 

By  virtue  of  optical  antenna  gain  reciprocity*,  the  ability  to 
calculate  the  temporal  statistics  of  the  variations  in  angle-of-arrival 
can  be  utilized  at  an  optical  receiver  aperture  to  determine  the  laser 
beam  pointing  Jitter  if  that  same  aperture  were  functioning  as  a laser 
transmitter  used  to  irradiate  what  before  was  the  optical  source.**  In 
this  analysis,  the  designator  will  be  considered  a receiver  and  the 
target  an  optical  source.  The  angle-of-arrival  temporal  statistics 
will  be  calculated  at  the  designator-" receiver". 

Fortunately,  a fairly  compact  set  of  theoretical  results  have 
recently  been  developed  for  the  subject  of  angle-of-arrival  temporal 


•Fried,  D.  L.  and  Yura,  H.  T,,  "Telescopic-Performance  Reciprocity 
for  Propagation  in  a Turbulent  Medium,"  J.  Opt.  Soc.  Am..  Volume  62. 
1972,  p.  600.  

••Hence,  H.  V.  and  Fried,  D.  L. , "experimental  Test  of  Optical 
Antenna-Gain  Reciprocity, 'J  J.  Opt,  Soc.  Am..  Volume  63,  1973,  p.  1015. 
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statistics.*  With  this  theory,  the  angle-of-arrival  statistics  can 
easily  be  formulated  for  a particularly  general  baseline  propagation 
problem.  It  is  then  a matter  of  geometric  argument  and  simple  algebraic 
manipulation  to  make  this  result  apply  to  the  propagation  geometry  of 
interest  in  the  case  of  the  laser  target  designator. 

In  the  following  section,  the  existing  theory  for  angle-of-arrival 
temporal  statistics  shall  be  reviewed.  In  the  next  section,  the  geometric 
arguments  will  be  provided  that  allow  translation  of  this  result  to  the 
geometry  of  the  laser  target  designator  problem.  After  that,  a sample 
numerical  problem  will  be  worked.  In  the  final  section  it  will  be  shown 
how  a suitably  designed  random  pointing  error  generator  might  be  config- 
ured in  digital  computer  simulation  to  produce  a signal  matching  the 
calculated  temporal  statistics. 


II.  PROPAGATION  FORMULATION 

A receiver  of  aperture  diameter  D,  receiving  laser  radiation 
of  wavelength  \,  is  considered  from  a point  source  at  a range  R.  The 
direction  from  the  point  source  to  (the  center  of)  the  receiver  is  taken 
to  be  the  z-axis,  letting  z ■ 0 at  the  point  source,  and  z * R at  the 
receiver  aperture.  This  geometry  is  depicted  In  Figure  1. 


APERTURE 

Figure  1.  Baseline  propagation  geometry. 


A receiver  aperture  of  diameter  D is  viewing  a point  source  at  range  R. 
The  z-axis  corresponds  to  the  direction  from  the  point  source  to  the 
aperture  (z  - 0 at  the  point  source).  For  the  baseline  problem,  cal- 
culation of  the  angle-of-arrival  power  spectrum  for  the  wavefront 
arriving  at  the  aperture  from  the  point  source  is  of  interest.  The 
effective  wind  is,  as  shown,  the  component  of  the  wind  (and  associated 
beam  slew  effects)  perpendicular  to  the  z-axis. 


^Greenwood,  D,  P.  and  Fried,  D.  L,,  Power  Spectra  Requirement  for 
Wavefront-Compensative  Systems.  Rome  Air  Development  Center,  AFSC,  Griff las 
Air  Force  Base,  New  York  13441,  September  1975,  RADC-TR-75-227. 


One  component  of  the  angle-of-arrival  variations  of  the  wavefront 
as  it  reaches  the  receiver  aperture,  measured  in  radians  shall  be  denoted 
by  a(t),  (Only  a single  component  is  necessary  because  the  statistics 
will  be  the  same  for  both  components  and  the  results  will  apply  equally 
well  to  each.)  It  can  be  shown  that  0£(t)  has  zero  mean  value  and  sta- 
tionary statistics.  Hence,  the  angle-of-arrival  temporal  coherence 
function  can  be  expressed  as 

C^Cr)  = <a(t)  a(t  + t)>  (1) 

where  t denotes  time,  measured  in  seconds.  The  associated  power  spec- 
trum can  be  written  as 

dr  cos  (2x  f t)  C (t)  (2) 

o 3 


F <f)  = 4 f 


with  f denoting  temporal  frequency,  measured  in  hertz.  The  fact  that 
C^Or)  is  an  even  function  insures  that  there  is  no  sine,  or  imaginary, 

part  to  the  power  spectrum.  The  calculation  of  F^Or)  will  be  the  basic 

concern  in  this  section. 


It  has  been  shown  by  Greenwood  and  Fried*  that  if  the  propagation 
range  z ■ 0 to  z - R is  subdivided  into  N adequately  small  regions,  of 
width  and  with  midpoints  at  z^,  where  i * 1 to  N,  then  the  following 

equations  can  be  written: 

F-<£>  ■ it  w 

where 

F (f)  - 1.32  x 10'2  (\/D) 2 (D/r  4)5/3f  f'2/3  G (f/f  ,). 

U,1  0,1  0,1  OC  0,1 

(4) 

The  function  G (f/n  .)  is  defined  by  the  expression 

Ur  V|1 

1 if  0 < f < 0.332  f , 

o,i 

G (f/ f .)  -{  1.12-0.361  (f/f  .)  if  0.332  f . < f < 3.10  f . 
a o,i  o,i  o,i  o,i 

0 if  3.10  f , < f (5) 

o,i 

*lbid. 


The  quanties  r . and  f . are  determined  by  the  local  turbulence  and 
°.t  o,i  th 

effective  wind  characteristics  in  the  i segment  of  the  propagation 

path.  With  C . ^ denoting  the  refractive-index  structure  constant  in 
th  N*1  th 

the  i segment,  V . denoting  the  effective  wind  velocity  in  the  i 
i ett,i 

segment,  z.  denoting  the  position  on  the  z-axis  of  the  midpoint  of  the 

t.  X 

i segment,  and  denoting  the  width  of  the  i segment,  the  following 
equations  can  be  written: 

ro,i  - l16-7  “l  C»,l2  («1/R)5/Vl'3/5  (6) 


and 


Veff.i 
n D (zi/R) 


(7) 


Equations  (3)  to  (7)  provide  a straightforward  basis  for  the  calculation 
of  angle-o f-arrival  fluctuations  for  a receiver  viewing  a point  source. 

This  consititutes  the  baseline  problem.  For  purposes  of  definite- 
ness, a computer  program  has  been  prepared  in  the  form  of  a subroutine 
for  the  calculation  of  the  baseline  power  spectrum.  This  is  listed  in 
Appendix  A. 


III.  LASER  TARGET  DESIGNATOR  PROPAGATION  GEOMETRY 

The  laser  target  designator  has  a rather  small  aperture 
diameter  compared  to  the  diameter  of  the  transmitted  laser  spot.  If 
the  laser  target  designator  diameter  is  and  the  angular  divergence 

of  the  laser  beam  is  9,  then  the  diameter  of  the  laser  spot  will  be 
Dg  - Dt  + 8 , (8) 

where  is  the  range  from  the  laser  target  designator  to  the  laser 

spot.  (The  Interest  here  is  restricted  to  the  gedanlten  situation  where 
the  laser  spot  is  defined  on  a virtual  screen  perpendicular  to  the 
z-axis.) 

To  calculate  the  angle  of  arrival  fluctuations  at  the  laser  spot, 
the  baseline  formalism  can  be  utilized,  except  that  it  would  be  assumed 
that  there  was  a point  source  located  a distance 
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H • V9 


(9) 


behind  the  laser  target  designator  aperture.  A further  assumption  is 
that  there  is  no  turbulence  in  this  added  region.  The  geometry  of  this 
configuration  is  depicted  in  Figure  2.  In  this  virtual  formulation, 
the  problem  posed  is  directly  reduced  to  a version  of  the  baseline  case. 


r 


TURBULENCE-FREE 

REGION 


A 


■\r 


TURBULENT  REGION 

A 


>1 


Figure  2.  Virtual  propagation  geometry  for  laser  target 
designator. 

The  laser  target  designator  projects  a laser  spot  at  range  which, 

because  of  the  laser  beam  spread,  is  considerably  larger  than  the  aper- 
ture diameter  of  the  laser  target  designator.  The  wavefront  angle-of- 
arrival  at  the  laser  spot  can  be  evaluated  by  treating  the  spot  as  a 
virtual  aperture  viewing  a virtual  point  source  located  behind  the  laser 
target  designator.  With  the  laser  target  designator  ignored,  and  treat- 
ing the  turbulence  as  being  restricted  to  the  R^  region  with  the  R^ 

region  being  turbulence- free,  the  problem  reduces  to  the  baseline  case. 


The  angle-of-arrival  can  be  analyzed  at  the  laser  spot  in  terms 
of  the  baseline  formulatio  by  doing  the  following: 


1)  The  baseline  parameter  D is  taken  to  be 

D - Dg  . (10) 

2)  The  baseline  parameter  R is  taken  to  be 

R - Rx  + R2  . (11) 
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3)  The  propagation  path  in  the  baseline  case  is  taken  to  be  sub- 
divided into  N + 1 regions  i = 0 to  N,  with  the  i = 0 segment  having 
a width 


Az^.  a R2  , for  i = 0, 


(12) 


and  a refractive  index  structure  constant 
„ 2 „ 


•'N.i 


for  i = 0 


(13) 


4)  For  the  part  of  the  propagation  path  corresponding  to  R , 

2 1 
the  values  of  N,  Az , and  C . are  taken  in  the  same  manner  if  the 
i N » 1 

propagation  path  were  just  over  the  length  R^. 


5)  The  values  of  z are  taken  to  be  those  which  would  be  obtained 
if  the  propagation  path  were  just  over  R^,  except  that  R2  is  added  to 
each  value;  thus. 


i-1 

z =* **  R„  + T Az  + 1/2  Az  . (14) 

j-1  J 

th 

6)  In  the  actual  wind  perpendicular  to  the  z-axis  in  the  i 
section  is  V.  , and  the  laser  target  designator  is  slewing  at  an 

Aii- 

angular  rate  0,  then  the  following  equation  results: 


ef  f ,i 


V , ± 0 


(zi  " 


R2> 


(15) 


where  the  plus  or  minus  sign  is  chosen  accordingly  as  the  angular  slew 

direction  is  opposed  to  or  parallel  to  the  sense  of  V 

A*1 


With  these  substitutions  and  using  the  algorithm  of  the  baseline 
case  as  represented  in  the  subroutine  of  Appendix  A,  the  power  spectrum 
of  each  component  of  the  angle-of-arrival  of  the  wave-front  can  be  cal- 
culated from  a point  source  as  it  reaches  the  laser  spot  position  in 
accordance  with  the  problem  definition  in  Figure  2.  it  is  a straight- 
forward matter  to  set  up  a computer  subroutine  for  this  calculation.. 
This  is  done  in  Appendix  B,  but  first  one  further  matter  must  be  noted. 


Having  determined  the  angle-of-arrival  fluctuation  statistics  at 

* 

the  laser  spot,  reciprocity  arguments  can  now  be  applied  the  same 
results  as  the  statistics  of  beam  pointing  variations  when  the  laser 
spot  is  a laser  transmitter  of  diameter  transmitting  a beam  focused 


*Op.  Cit.,  Fried  and  Yura. 

**0p.  Cit.,  Hence  and  Fried 
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on  the  position  of  the  virtual  point  source  at  range  + R^.  The  power 

spectrum  of  pointing  jitter  at  the  laser  spot  position  is  identical 
with  what  would  be  calculated  for  the  angle-of-arrival  spectrum.  Now, 
considering  the  laser  spot  position  as  the  transmitter,  the  angle-of- 
arrival  spectrum  of  the  partially  converged  but  not  fully  focused  beam 
can  be  investigated  when  it  reaches  the  position  of  the  laser  target 
designator.  It  can  be  shown  through  any  of  several  arguments  that  the 
angle-of-arrival  power  spectrum  at  the  laser  target  designator  position 
2 

is  just  (D^/D^)  times  the  beam  pointing  jitter  and  thus  times  the 

angle-of-arrival  jitter  at  the  laser  spot  position.  The  simplest  argu- 
ment for  this  follows  from  the  LaGrange  invariant  which  requires  that 
in  a passive  optical  system  (which  is  what  the  turbulent  atmosphere  is 
at  any  instant)  angle  times  aperture  (i.e.,  beam)  diameter  at  one 
position  in  the  system  must  equal  angle  times  aperture  (i.e.,  beam) 
diameter  at  any  other  position  in  the  system.  Because  the  angle-of- 
arrival  power  spectrum  is  proportional  to  angle  squared,  the  ratio  of 

2 

diameters  squared  is  utilized;  hence  the  (Dg/D^)  factor. 

It  is  immediately  obvious  from  reciprocity  considerations  that  the 
pointing  jitter  induced  turbulence  in  the  laser  target  designator  trans- 

2 

mitting  a beam  to  the  laser  spot  is  also  equal  to  (DS^DT)  times  the 

angle-of-arrival  spectrum  at  the  laser  spot  position  of  the  wavefront 
from  the  virtual  point  source:  The  formalism  for  calculating  this  latter 

quantity  was  set  up  in  the  first  part  of  this  section.  In  Appendix  B, 
this  formalism  is  utilized  to  allow  calculation  of  the  laser  target 
designator  turbulence  induced  pointing  jitter  power  spectrum. 


IV.  SAMPLE  RESULTS 


With  the  capabilities  defined  by  the  programs  in  Appendices  A 
and  B,  a set  of  sample  results  can  be  calculated.  In  this  section,  first 
a sample  calculation  of  the  power  spectrum  shall  be  performed  and  then 
it  will  be  shown  how  this  results  can  be  used  to  generate  an  appropriate 
random  sample  of  a time  varying  laser  target  designator  pointing  jitter. 


A problem  is  considered  in  which  the  laser  target  designator  is 

operating  at  a wavelength  of  \ = 1.06  X 10  ^ m,  has  an  aperture  diameter 

3 

of  D ■ 0.04  m,  and  is  illuminating  a target  at  a range  R = 3 x 10  m. 

T _4  1 

The  laser  beam  spread  (full  angle)  is  9 = 5 X 10  rad;  the  target  is 
assumed  to  be  moving  at  a crossing  velocity,  i.e.,  a velocity  perpen- 
dicular to  the  line  of  sight  of  25  km/hour  (6.94  m/sec).  At  a 3-km 
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• -3 

range,  this  gives  rise  to  a beam  slewing  angular  rate  of  0 = 2,31  X 10 
rad/sec.  It  is  assumed  that  the  propagation  is  between  two  hills  and 
across  a valley  so  that  the  strength  of  turbulence  is  lowest  at  the 
center  of  the  path  and  approximately  equally  higher  at  both  ends;  the 
wind  velocity  is  lowest  at  the  two  ends  and  highest  at  the  center. 


The  propagation  path  shall  be  considered  to  be  subdivided  into 
N = 10  equal  length  segments,  so  that  the  segment  length  are  Az^  = 300  m 

(i  = 1,2,,..,  10).  The  refractive-index  structure  constant  is  taken  to 

-14  -2/3  -14 

vary  between  10  X 10  m at  the  two  ends  of  the  path  to  2 X 10 

-2/3 

m at  the  center  of  the  path.  Explicitly,  the  following  equations 
are  set: 


10  x 

7 X 
5 X 
3 X 
2 X 


10 

10 

10 

10 

10 


-14 

-14 

-14 

-14 

-14 


» 


» 


» 


» 


(16) 


The  actual  wind  velocity  component  perpendicular  to  the  propagation 
path  is  taken  to  vary  between  10  m/sec  near  the  center  of  the  path  to 
3 m/ sec  near  the  two  ends  of  the  propagation  path  (l.e.,  between  approx- 
imately 22  mph  and  7 mph).  Explicitly,  the  following  equations  are  set: 


Vx.i-Vx,io 
VX,2  “ VX,9 
VX,3  " VX,8 


3 m / sec  , 

6 m/sec  , 

8 m/sec  , 

9 m/  sec  , 

10  m/ sec  . 


These  values  completely  define  the  propagation  problem. 


(17) 


A computer  run  has  been  set  up  using  the  DATA  instructions  shown 
in  Appendix  C to  calculate  the  laser  target  designator  atmospheric 
turbulence  induced  pointing  jitter  for  the  M ■ 9 (logarithmically 
uniformly  spaced)  frequencies  of  f - 0.1,  0.1778,  0,3162,  0.5623,  1,000, 
1,778,  3.162,  5.623,  and  10.00  Hz.  The  results  of  this  calculation 
are  shown  in  Appendix  C and  graphed  in  Figure  3. 
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POINTING  JITTER  POWER  SPECTRUM.  F „(f)  (r«T/Hi) 


Figure  3.  Later  target  designator  pointing  jitter  power  spectrum. 

Pointing  jitter  is  induced  by  the  turbulence  and  wind  speed  distri- 
bution given  as  an  example  in  the  test.  The  complete  set  of  problem 
parameters  are  specified  by  the  "DATA"  statements  given  in  Appendix  C. 
The  computation  of  the  power  spectrum  was  carried  out  using  the  pro- 
grams in  Appendices  A and  B. 
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V.  RANDOM  JITTER  GENERATION 


The  problem  now  is  to  find  a means  of  easily  generating  a 
random  sequence  with  this  power  spectrum.  It  is  suggested  that  the 
simplest  way  to  do  this  is  to  define  a set  of  simple  digital  R-C  low 
pass  6-dB/octave  filters  which  when  summed  will  have  a transfer  function 
equal  to  the  square  root  of  the  pointing  jitter  power  spectrum.  If  a 
random  white  noise  is  put  into  such  a network,  the  output  will  be  random 
noise  with  the  appropriate  power  spectrum.  The  matter  of  scaling  will 
be  discussed  shortly,  but  first  the  generation  in  a computer  simulation 
of  the  random  white  noise  and  the  formulation  of  the  appropriate  R-C 
network  in  the  computer  will  be  discussed. 

For  all  of  the  computer  work,  a At  ■ 0.001  sec  time  increment  will 

be  utilized.  This  means  that  frequencies  up  to  (2At)  ^ = 500  Hz  can  be 
handled  and  frequencies  up  to  approximately  100  Hz  can  be  handled  ade- 
quately.* A random  sequence  of  gaussianly  distributed  independent  random 
values  with  unity  variance  is  generated  associating  each  one  with  a 
sequence  of  time  intervals  At,  then  the  resultant  signal  will  appear  to 
be  white  noise  with  unity  variance  and  a power  spectrum  that  is  abruptly 

truncated  at  the  Nyquist  frequency,  f^  *=  (2At)~*.  Because  the  variance 

of  a random  signal  equals  the  integral  of  the  power  spectrum  over  the 
frequency  range  0 to  °°,  it  follows  that  the  power  spectral  density  of 
this  white  noise  is 

N “ 1/f  - 2At  . (18) 

o N 

A single  R-C  low  pass  digital  filter  is  digitally  implemented  in 

accordance  with  Figure  4.  Here  the  box  labeled  "z  denotes  a time 
delay  between  input  and  output  of  on«  cycle  time,  At.  If  the  factor  k 
is  chosen  to  be 

k - exp  (-2x  fR(,  At)  , (19) 

and  the  factor  K is  chosen  to  be 

K - 1 -k  , (20) 

then  in  the  limit  of  sample  data  approximations,  the  output  signal  will 
appear  to  be  the  input  signal  filter  by  an  R-C  low  pass  filter  with  the 
knee  of  its  transfer  function  at  f__,  It  is  noted  that  this  filter  has 


*If  higher  frequencies  are  required,  it  is  merely  necessary  to  use 
a smaller  value  for  At. 
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unity  gain,  i.e.,  if  a dc  signal  is  input,  the  output  will  be  a dc  signal 
of  the  same  amplitude.  The  transfer  function  of  this  filter  will  be 
denoted  by  Y(fon;f).  A set  of  these  filters  will  be  summed;  each  with 

separate  f_r  and  separate  weighting  factor,  A,  so  as  to  obtain  a com- 
posite  transfer  function  Y,p  Qt  (f),  whose  magnitude  squared  will  approx- 
imate the  pointing  jitter  power  spectrum  in  Figure  4. 


Figure  4.  Recursive  filter  design  for  digital  computer  simulation 
of  an  R-C  low  pass  filter. 

The  box  labeled  z denotes  a single  time  delay  of  duration  At,  To 
obtain  a unity  gain  filter  with  the  knee  of  its  transfer  function  at 
fRC’  the  multiPlicative  constants  k and  K are  set  equal  to  exp  (-2rt 

fRC  At)  and  1 - exp  (-2jt  fR(,  At)  , respectively. 

Thus,  the  following  equation  can  be  written: 

J 

^T  't  ^ ^ » (21) 

«J®  1 

and 

l*Tot  * Power  Spectrum  of  Figure  3 . (22) 

The  value  of  J is  chosen  to  allow  reasonably  accurate  approve  mat ion  in 
Equation  (22),  (Care  is  required  in  the  choice  of  the  A/s  and  L 's 
to  insure  a good  approximation.)  J ~J*J 

Now  if  the  sequence  of  unity  variance  random  values  is  multiplied 

by  Nq  and  passed  through  the  filter  YIot  (f),  the  output  will  have 

the  statistics  of  the  laser  target  designator  turbulence  induced  point- 
ing jitter. 
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Appendix  A.  CALCULATION  OF  THE  BASELINE  POWER  SPECTRUM 


1000 

1001- 

1002 

1003 

1004 

1005 

1006 

1007 

1008 

1009 

1010 
101 1 
1012 

1013 

1014 

1015 

1016 

1017 

1018 

1019 

1020 
1021 
1022 

1023 

1024 
1050 
1055 
1060 
1065 
1070 
|07S 
1080 
1085 
1090 
109S 


REM 

REM  THIS  P Rd  (.  RAM  . W R I T EN  IN  BASIC  WILL  CALCULATE 
REM  THE  POWER  SPECTRUM  OF  THE  ANC LE-0 F - ARR I VAL  FOR 
REM  THE  BASELINE  CASE  OF  AN  APERTURE  VIEWING  A POINT 
REM  SOURCE.  THE  PROGRAM  FUNCTIONS  AS  A SUBROUTINE. 

REM  AND  IS  ENTERED  V.  ITH  THE  FOLLOWING  PARAMETERS  SETl 
REM 

REM  LsLASER  WAVELENGTH.  LAMBDA 

REM  D*APERTURE  DIAMETER 

REM  R=PRO  PA GAT  ION  PATH  LENGTH 

REM  F=F REQUENCY  FOR  WHICH  THE  POWER  SPECTRUM 

REM  IS  T0  RE  EVALUATED 

REM  N=NUMBER  OF  SEGMENTS  IN  THE  PROPAGATION 

REM  PATH 

REM  Z t f ) = RANGE  FROM  BASELINE  CASE  POINT  SOURCE 

REM  TO  THE  MIDPOINT  OF  THE  I-TH  SEGMENT 

REM  Zl C I ) “LENGTH  OF  THE  I-TH  SEGMENT 

REM  C<I)=C-SUB-N  SQUARED  FOR  THE  I-TH  SEGMENT 

REM  V< I )*V-SUB-EFF  FOR  THE  I-TH  SEGMENT 

REM 

REM  THE  SUBROUTINE  IS  EXITED  WITH 
REM 

REM  FI =VALUE  3 F THE  POWER  SPECTRUM  FOR  FREQUENCY  F. 

REM 

REM 

LET  F 1 =0 
FOR  1=1  TO  N 

LET  FO=VC I J/C3.  1 41 59*D*ZC II/R) 

LET  R0=< 16.7*Z1 t I J*Cl  Il*<Zm/R)»  1 . 66667/ CL*L>  >»  (-.6) 

IF  F«= • 332*F0  LET  G=! 

IF  F».332*F0  IF  F<=3. 1*F0  LET  G= I . 1 2- . 36 1 *F/TO 
IF  F>3. 1*F0  LET  6=0 

LET  FI *F1 ♦ 1 . 32E-02*<L/D>* 2*<  CD/RO)»  5/<F*F*F0> >t . 333333*G 

NEXT  I 

RETURN 
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Appendix  B.  CALCULATION  OF  THE  POWER  SPECTRUM  OF  EACH  COMPONENT 


100  REM 

101  REM  THIS  PROGRAM, WRITER  IN  BASIC  WILL  CALCULATE  THE 

102  REM  POWER  SPECTRUM  OF  EACH  COMPONENT  OF  ATMOSPHERIC 

103  REM  TURBULENCE  INDUCED  POINTING  JITTER  OF  A LASER 

104  REM  TARGET  DESIGNATOR.  THE  PROGRAM  MAKES  USE  Or  THE 

105  REM  SUBROUTINE  IN  APPENDIX  A.  THE  PROGRAM  IS  STARTED 

106  REM  WITH  THE  PARAMETERS  DEFINING  THE  PROBLEM  SITUATION 

107  REM  AND  THE  NATURE  OF  THE  REOU I RED  RESULTS  GIVEN  BY 

108  REM  A SET  OF  ••DATA"  STATEMENTS.  AS  LISTED  HERE  THESE 

109  REM  "DATA"  STATEMENTS  ARE  FILLED  WITH  ZEROS.  FOR  ACTUAL 

110  REM  PROGRAM  OPERATION  THESE  ZERO  VALUES  HAVE  TO  BE 

111  REM  REPLACED  BY  THE  APROPRI ATE  VALUES. 

112  REM 

113  REM  LINE  150  IS  A "DATA"  STATEMENT  LISTING  IN  ORDER 

114  REM  THE  VALUES  OF  THE  FOLLOWING! 

115  REM  L=LASER  WAVELENGTH,  LAMBDA  (IN  METERS) 

116  REM  D1=LASER  TARGET  DESIGNATOR  APERTURE  DIAMETER, 

117  REM  D-SUB-T  (IN  METERS) 

118  REM  T»LASER  BEAM  SPREAD,  THETA  (IN  RADIANS) 

119  REM  TJaLASER  BEAM  SLEW  RATE,  THETA-DOT  (IN 

120  REM  RADIANS  PER  SECOND) 

121  REM  R1*PR0PAGATI0N  RANGE  FROM  THE  LASER  TARGET 

122  REM  DESIGNATOR  TO  THE  LASER  SPOT,  R-SUB-1 

123  REM  < IN  METERS) 

124  REM 

ISO  DATA  0,0, 0,0,0 

155  READ  L 

156  READ  D! 

157  READ  T 

158  READ  T! 

159  READ  R1 

160  LET  D2=D1*T*RI 

161  LET  R2=Dl/7 

165  REM 

166  REM 

167  REM 

168  REM 

169  HEM 

170  REM 

171  REM 

172  REM 

173  REM 

174  REM 

175  REM 

176  REM 

177  REM' 

178  REM 

179  REM 

180  REM 

181  REM 

182  REM 


LINE  190  AND  THE  IMMEDIATELY  FOLLOWING  LINES  ARE 
"DATA”  STATEMENTS  GIVING  IN  ORDER  THE  VALUES  OF 
THE  FOLLOWING) 

N*NUMBER  OF  SEGMENTS  THE  PROPAGATION  PATH  IS 
CONSIDERED  TO  BE  DIVIDED  INTO  FOR  THE 
COMPUTATION 

ZKI)aA  SET  OF  N-VALUES  CORRESPONDING  TO  THE 
LENGTH  OF  THE  I-TH  SEGMENT,  DELTA-Z- 
SUB-I  CIN  METERS) 

C(I)«A  SET  0F  N-VALUES  OF  THE  REFRACTIVE- INDEX 
STRUCTURE  CONSTANT  IN  THE  I-TH  SEGMENT, 
C-SUB-N-SOUARED-SUB- I (IN  METERS* 2/3) 

V I ( I }»A  SET  OF  N-VALUES  OF  THE  ACTUAL  WIND  VELOCITY 
PERPENDICULAR  TO  THE  PROPAGATION  DIRECTION 
IN  THE  I-TH  SEGMENT,  V-SUB-PERPENDICULAR 
(IN  METERS  PER  SECOND) 
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190  DATA  0 

19!  DATA  0#0#0»0*0*0*0*0»0*0 

192  DATA  0#0*0«0#0*0«0«0#0#0 

193  DATA  0#0#0«0#0*0«0#0#0#0 
200  READ  N 

205  DIN  Z1CN]#ZCNI*CCN]#V£NI 

210  LET  Z=R2 

215  LET  Zl*0 

220  FOR  1*1  TO  N 

225  READ  Z1CI3 

230  LET  Z*Z«-.S*(Z1*Z1CI3> 

235  LET  ZCIJ*Z 
240  LET  Z1*Z1CIJ 
245  NEXT  I 
250  FOR  1=1  TO  N 
255  READ  CCI3 
260  NEXT  I 
265  FOR  1*1  TO  N 
270  READ  VI 

275  LET  VCIJ«V1<*T1*(Z(IJ-R2: 

200  NEXT  I 

285  REM 

286  REM  LINE  300  AND  THE  LINES  IMMEDIATELY  FOLLOWING  ARE 

287  REM  "DATA**  STATEMENTS  SPESIFYING  THE  FOLLOWING/ 

288  REM  M*THE  NUMBER  OF  FREQUENCIES  FOR  WHICH  THE 

289  REM  POINTING  JITTER  POWER  SPECTRUM  IS  TO 

290  REM  CALCULATED 

291  REM  F C I )*A  SET  OF  M-VALUES  CORRESPONDING  TO  THE 

292  REM  VARIOUS  FREQUENCIES  FOR  WHICH  THE  POWER 

293  REM  SPECTRUM  IS  TO  BE  CALCULATED  (IN  HERTZ>* 

294  REM 

300  DATA  O 

301  DATA  0«0*0#0#0«0#0#  0*0*0 
305  READ  M 

310  DIM  FCMI 
315  FOR  1*1  TO  M 
320  READ  FLU 
325  NEXT  I 

330  REM 

331  REM  ALL  OF  THE  PROBLEM  PARAMETERS  ARE  ESTABLISHED  AT 

332  REM  THIS  POINT*  THE  PROGRAM  NOW  SETS  UP  THE  VARIOUS 

333  REN  PARAMETERS  TO  CALL  THE  SUBROUTINE#  AND  AFTER  EACH 

334  REM  CALL  PRINTS  OUT  THE  RELEVANT  RESULTS* 

335  REM 

340  LET  D*D2 
345  LET  R*R1*R2 
350  FOR  J*1  TO  M 
355  LET  F*FCJ) 

360  GO SUB  1000 

365  LET  Fl*Fl*(D2/Dl >»2 

370  PRINT  "FREQUENCY  • "iFl  TAB  (25>I"P0WER  SPECTRUM  ■ "IF I 
375  NEXT  J 
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Appendix  C.  CALCULATION  OF  THE  LASER  TARGET  DESIGNATOR 
ATMOSPHERIC  TURBULENCE 


In  Che  text,  a sample  problem  was  set  up  for  evaluation.  Here, 
the  appropriate  DATA  statements  for  the  problem  and  then  the  computed 
power  spectrum  are  exhibited  in  place  of  the  empty  DATA  statements  of 
Appendix  B as  follows: 

150  DATA  1.06E-6*4E-2*5E-4*2.31E-3*3E3 


190  DATA  10 

191  DATA  300#  300# 300# 300# 300#  300* 300*  300#  300#  300 

192  DATA  1E-13.7E-14* 5E-14*3E-14* 2E-14 

193  DATA  2E-|4*3E-14*  5E-14*7E-14#  IE- 13 

194  DATA  3*6, 8*9#  10*  i0#9#8#6#3 


300  DATA  9 

301  DATA  0*1*0. 1778*0. 3162*0. 5623*  1. 0*1. 778*3. 162#S. 623*10*0 

Line  150  gives  \,DT>  9,9,  and  , in  that  order.  Line  190  says  that 

the  propagation  range  is  divided  into  10  intervals.  Line  191  gives  the 
value  of  AZj  for  each  of  the  10  intervals.  Line  192  gives  the  value  of 

^or  ®ach  of  the  10  intervals.  Line  193  gives  the  values  of  ^ 

for  each  of  the  10  intervals.  Line  300  says  that  the  power  spectra  will 
be  computed  for  9 frequencies.  Line  301  gives  the  values  of  the  9 
frequencies . 

Using  the  programs  in  Appendices  A and,  B,  and  the  preceding  DATA  state- 
ments, the  following  results  were  calculated.  These  results  are  plotted 
in  Figure  3. 
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